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The influence of non-linear interaction effects on the characteristics of hadronic 
collisions and on the development of extensive air showers (EAS) is investigated. 
Hadronic interactions are treated phenomenologically in the framework of Gri- 
bov's Reggeon approach with non-linear corrections being described by enhanced 
(Pomeron-Pomeron interaction) diagrams. A re-summation algorithm for higher 
order enhanced graphs is proposed. The approach is applied to develop a new 
hadronic interaction model QGSJET-II, treating non-linear effects explicitely in 
individual hadronic and nuclear collisions. The model is applied to EAS modeliza- 
tion and the obtained results are compared to the original QGSJET predictions. 
Possible consequences for EAS data interpretation are discussed. 



1. Introduction 

Currently the physics of hadronic interactions remains one of the most in- 
triguing fields both for experimental and for theoretical research. Being 
the subject of investigations in collider experiments, it has a different role 
in high energy cosmic ray (CR) studies. There, a proper understanding of 
hadronic interaction mechanisms is of vital importance for a correct inter- 
pretation of CR data, for a reconstruction of energy spectra and of particle 
composition of the primary cosmic radiation, and finally, for inferring in- 
formations on the sources of CR particles, whose energies extend to ZeV 
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region, exceeding by far those attainable at man-made machines. 

Despite a significant progress in QCD over the past two decades, only 
phenomenological treatment is generally possible for minimum-bias hadron- 
hadron (hadron-nucleus, nucleus-nucleus) collisions. In particular, Gribov's 
Reggeon scheme 1 proved to be a suitable framework for developing success- 
ful model approaches, like Quark-Gluon String or Dual Parton models, 2 
which in turn provided a basis for a number of Monte Carlo (MC) gener- 
ators of hadronic and nuclear collisions, extensively used in accelerator or 



Still, all the mentioned MC models are based on a linear picture, the 
elastic scattering amplitude being defined by (quasi-)eikonal contributions 
of independent Pomeron exchanges. Meanwhile, in the limit of very high 
energies one expects a significant increase of parton density in colliding 
hadrons (nuclei), which gives rise to essential non-linear interaction effects. 7 
On the other hand, in Gribov's scheme the latter are traditionally described 
as Pomeron-Pomeron interactions. 8,9 Here we propose a re-summation pro- 
cedure both for enhanced contributions to the elastic scattering amplitude 
and for various unitarity cuts of the corresponding diagrams. This gives us 
a possibility to account for non-linear corrections when calculating hadron- 
hadron cross sections and to develop a new MC model, QGSJET-II, 10 which 
treats non-linear effects explicitely in individual hadronic and nuclear col- 
lisions. The model is applied to simulate extensive air showers (EAS), 
induced by CR particles in the atmosphere, which allows to investigate the 
influence of non-linear corrections on the air shower characteristics and to 
draw possible consequences for EAS data interpretation. 

2. Linear Scheme 

Gribov's Reggeon approach 1 describes a high energy hadron-hadron colli- 
sion as a multiple scattering process, where elementary re-scatterings are 
treated phcnomcnologically as Pomeron exchanges - Fig. 1. Thus, hadron a 
- hadron d elastic amplitude can be obtained summing over multi-Pomeron 
exchange contributions' 12,11 (see also Ref. 12): 



CR fields, e.g., DPMJET, 3 VENUS, 4 or QGSJET. 5,6 




dx+dx^\«\ d k G r ad (x+xis,b) 



1 > ■••> x n )i 



(1) 



a In the high energy limit all amplitudes can be considered as pure imaginary. 
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Figure 1. A general multi-Pomeron contribution to hadron-hadron scattering ampli- 
tude. Elementary scattering processes (vertical thick lines) are described as Pomeron 
exchanges; thin lines correspond to constituent partons, to which Pomerons are coupled. 



where s and b are the cm. energy squared and the impact parameter 
of the interaction, G r ad (x + x~s, b) is the unintegrated Pomeron exchange 
eikonal (for fixed values of Pomeron light cone momentum shares x+ ) , and 
, x n ) is the light cone momentum distribution of constituent par- 
tons - "Pomeron ends" (here, quark-antiquark pairs). Cj, A" are the rel- 
ative weights and strengths of diffractive eigenstates of hadron a in the 
Good- Walker formalism; 13 T,j C j = h ^jCjtf = l - In particular, a 
two-component model (J — 1,2) with one passive component, Af = 0, cor- 
responds to the quasi-eikonal approach, where A" = 1/Cf is the shower 
enhancement coefficient. 2 

Eq. (1) can be greatly simplified assuming that the integral over light 
cone momenta can be factorized: b 



J]_dx+dx l G v ad {xjx l s,b) 



i=i 



Njr\xt t ...,x+)N^(xl,...,x-) 



J dx+dx- G p ad (x + x- S ,b)Nj 1 1 \x + )N d :i \x-) 



This leads to traditional eikonal formulas: 

f ad (s,b)=iJ2 C J C £ [l-e-^^ s V 



X r ad (s,b) = / dx+dx- G p ad (x + x-s,b) iV«(x + ) N™(x-) 



(2) 



(3) 



(4) 



In this scheme the Pomeron is an effective description of a microscopic 
parton cascade which mediates the interaction between the projectile and 
the target hadrons. It is convenient to divide the latter into two parts: "soft" 



3 Here we neglect energy-momentum correlations between multiple re-scatterings. 14 
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cascade of partons of small virtualitics \q 2 \ < Q 2 , and a perturbative parton 
evolution at \q 2 \ > Qq, where Qq is some cutoff for pQCD being applicable. 
Correspondingly, a "general Pomeron" will consist of two contributions: 
"soft" Pomeron for a pure non-perturbative process (all \q 2 \ < Qq) and a 
"semi-hard Pomeron" for a cascade which at least partly develops in the 



high virtuality region (some \q | > Qq) 



2\6,12,15 



Fig. 2: 



(5) 



The advantage of such a procedure is that in very high energy limit the 
"general Pomeron" is dominated by its semi-hard component and the en- 
ergy dependence of the eikonal G^ d (s, b) is governed asymptotically by the 
perturbative QCD evolution. 





soft Pomeron 



QCD ladder 



soft Pomeron 



Figure 2. A "general Pomeron" (l.h.s.) consists of the "soft" and the "semi-hard" 
Pomerons — correspondingly the 1st and the 2nd contributions on the r.h.s. 



The soft Pomeron eikonal can chosen in the usual form: 2 



la Id (s/ S0 ) a 'W- 

ap(0) ln- 



m 



—e 



4 ( i? c 



+ Ri+c'{0) In 



(6) 



where sq ~ 1 GeV 2 is the hadronic mass scale, ap(0) and ap(0) are the 
intercept and the slope of the Pomeron Regge trajectory, and j a , Ra are 
the coupling and the slope of Pomeron-hadron a interaction vertex. 

The dominant contribution to the semi-hard Pomeron comes from hard 
scattering of gluons and sea quarks c and can be represented by a piece 
of QCD ladder sandwiched between two soft Pomerons d - Fig. 2; the 



c For brevity, hard scattering of valence quarks is not discussed explicitely. 

d In general, one may consider an arbitrary number of ^channel iterations of soft and 

hard Pomerons. 16 
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corresponding eikonal may be written as 6,12. 




xG F a r( So /x+,b') G r if( S «/x- hl \b- b'\) 



(7) 



Here afj (ie^x^s, Q 2 ,) is the contribution of the parton ladder with the 
virtuality cutoff \q 2 \ > Qq; I, J and x~t,x7 are the types and the rela- 
tive light cone momentum shares of the ladder leg partons. The eikonal 
G'y ft (s, 6'), describing parton / momentum and impact parameter distri- 
bution in the soft Pomeron at virtuality scale Qq, is defined by Eq. (6), 
neglecting the small slope of Pomeron-parton coupling R 2 ~ 1/Qq ~ 
and using a parameterized Pomeron-parton I vertex ji(xh), Xh — So/s. 
Usual hadronic parton momentum distribution functions (PDFs) can be 
obtained convoluting G*f ft with the constituent parton distribution 
and integrating over the parton impact parameter b': 



Thus, the contribution of semi-hard processes to the integrated eikonal 
(4), integrated over the impact parameter b, can be written 12 ' 15 as a con- 
volution of hadronic PDFs fi/ a u/d){ x i Q 2 ) (obtained by evolving the input 
PDFs (8) from Qq to Q 2 ) with the parton scatter cross section dafp 2 /dp 2 : 



with pt being parton transverse momentum in the hard process, Mp - the 
factorization scale (here Mp — p 2 /4), and with the factor K accounting 
for higher order corrections. Clearly, the integrand in the curly brackets on 
the r.h.s. of Eq. (9) defines inclusive jet production cross section. 

Knowing the elastic scattering amplitude, Eq. (3), one can calculate 
total and elastic cross sections, elastic scattering slope for the interaction, 
etc. Moreover, cutting the corresponding diagrams of Fig. 1 according 
to the Abramovskii-Gribov-Kancheli (AGK) cutting rules 17 and collecting 
contributions of cuts of certain topologies, one can obtain relative weights 
for various configurations of the interaction. 2 The latter gives a possibility 
to develop MC generation procedures for hadron-hadron, hadron-nucleus, 
and nucleus- nucleus collisions. 5,6 ' 15 




(8) 




^fi/a(xt,M 2 F ) f J/d (x-,M 2 ) Q(M 2 - Ql)} , (9) 
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3. Pomeron-Pomeron Interactions 

The above-described linear picture cease to be valid in the "dense" regime, 
i.e. in the limit of high energies and small impact parameters of the inter- 
action. There, a large number of elementary scattering processes occurs 
and corresponding underlying parton cascades largely overlap and inter- 
act with each other. Such effects are traditionally described by enhanced 
diagrams. 8 ' 9 To develop a dynamic scheme it is assumed that Pomeron- 
Pomeron interactions are dominated by partonic processes at compara- 
tively low virtualities, 10,12 \q 2 \ < Qq, and that Pomeron-Pomeron coupling 
at larger virtualities becomes important only after reaching the satura- 
tion regime for parton densities at the scale Q\. Thus, we develop a 
scheme where multi-Pomeron vertexes involve only interactions between 
soft Pomerons or between "soft ends" of semi-hard Pomerons - Fig. 3. 




Figure 3. Contributions to the triple-Pomeron vertex from interactions between soft 
and semi-hard Pomerons. 

Basically we shall stay close to the 7r-meson dominance picture, 9 where 
all the vertexes g mn for the transition of to into n Pomerons have been 
expressed via a single additional constant r^p and where an asymptotic 
re-summation procedure has been proposed. 18 In particular, for the lowest 
in contribution with only one multi-Pomeron vertex (Fig. 4) one can 
obtain, summing over to, n > 1 and subtracting the term with m = n = 1 
(Pomeron self-coupling): 18 

-A? \ d k xlAso/xr, b 1 ) xUxrs, \b- b'\)} , (10) 
with Xap(*'^') = Xa7r(sjfr')j 7p = ln\ here and below the indexes j,k refer 
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Figure 4. Lowest order enhanced graphs; constituent partons not shown explicitely. 

to the diffractive eigenstates of hadrons a, d correspondingly. It is easy to 
see that in the large s, small b limit &Xad : ( s > ^) * s dominated by the 
last term in the integrand, i.e. by the subtracted self-coupling contribution. 
Therefore, asymptotically it was sufficient to consider a small sub-set of 
enhanced graphs, which can be obtained from the one in Fig. 4 iterating 
multi-Pomeron vertex in i-channel. In the "dense" limit re-summation of 
those diagrams reduces to the sum over subtracted Pomeron self-couplings 
and, after neglecting the slope of the multi-Pomeron vertex, Rp ~ 0, results 
in a re-normalization of the Pomeron intercept, 18 ap(0) = aqp(0)— Air r^hv- 
In this work we also assume an eikonal structure of the multi-Pomeron 
vertexes, 7p™p™^ = r 3 r 7™ +n ~ 3 , but treat 7r as a free parameter and neglect 
the vertex slope, Rp ~ 0. On the other hand, we neglect a momentum 
spread of "Pomeron ends" in the vertexes and define the eikonal x'p as 

xUS,b')= ( dx+ Ni 1 \x+)GUx+S 1 b'), (11) 

J so/s 

where G r aV is given by Eqs. (5-7), with d — > P. Correspondingly, for a 
Pomeron exchanged between two internal vertexes, e.g., in Fig. 5, we just 
use Gp P (so xp k /xp k+1 , \bk — &fc+i|), the latter being defined by Eqs. (5-7) 
with ad -> PP. 




Figure 5. An example of a "loop" graph. 
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Our goal is to develop a re-summation procedure which assures a smooth 
transition between the "dilute" (small s, large b) and "dense" limits and 
accounts for all essential contributions in both cases. Thus, we shall only 
neglect "loop" diagrams (Fig. 5) and "chess-board" graphs with more than 
three vertexes in a horizontal row - Fig. 6. The neglected diagrams give 




Figure 6. An example of a "chess-board" graph; only three vertexes are shown ex- 
plicitely. 



small contribution at low energies and/or large impact parameters, being 
proportional to large powers (> 2) of r^p. On the other hand, in the "dense" 
limit they are strongly suppressed by exponential factors. 18 For the example 
of Fig. 5, summing over any number ri\ of Pomerons exchanged between 
the target and the upper multi-Pomeron vertex (for the other vertex the 
procedure is similar), we can obtain a factor 



E 

711=0 



-^fcXdP^PiS.lb-bil) 



mi 



(12) 



We start with obtaining a "fan" diagram contribution Xa ^ via a re 
cursive equation - Fig. 7: 




>=2. 



+ 



x, b' 




x, b' 



Figure 7. Recursive equation for a "fan" diagram contribution Xa \S,b'), s = sq/x. 
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x f r ij \s,b')=^ aJ ,(s,b') 



T3P f 1 dXW 

X- 7r J s -f. xri 



i <FhG%(x Pl s,\b' -h\ 



E 



+ 1^3 I ^ G$*i**J, \b' ~ &il) 

A°7^ J 20. x Vl J 

1 _ e -\^»»(aoK,« _ A^xf^W^r A)l (13) 



Also we introduce a "generalized fan" contribution x at ^ 
recursive equation - Fig. 8: 



G-fan(j,fc) 



G-fan(j,fc)/ s T*/ 



(s,fe', s ,6) = x ys,6') + - 



^3P 



1 C?Xf 



\dH x 



(s /x Pl ,6i,s,6)| 



X _ e - A ? X a<J fa °"' f! '(so/xp 1 ,6i,s,6) 



(14) 



The difference between Xad fa fe ' and Xo™ (Fig- 9) is due to vertexes 
with both "fans" connected to the projectile and ones connected to the 
target. 




x,b' 




Figure 8. Recursive equation for a "generalized fan" x?j tan ^' k \s,b' , s,b), s = sq/x. 



Now we can obtain total eikonal for an "elementary scattering" process, 
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Figure 9. Symbolic notation for the difference between the "generalized fan" and "fan" 
contributions, \ ad — Xa ■ 



including all essential enhanced diagram contributions 10 - Fig. 10: 



r3V 



A?A^ 7 | 



fd 2 b' \(l 

o X V J l\ 



- p- X i *ad 



(l - e -X ) _ X a xtx O-^ x g-*»W) 



fan(fcj) fan(fc) 
Xd 



1 _ e - X i Xad 



2 (^"^ a 

-*tx ! r W \o G-fB.n(j,k)\ K ( G-fan(fc,j) P \ 
~ A J 'Xad J + -y ^Xda ~ XdV J 

U) -^Xa nU) )} - [a,j»d,k]}, (15) 



1 - e~ A ? x » 



where x Q rf — X a d {So/Xr, o , s, b), Xa — Xa (s /x P ,b ), 



Xap = Xap( s o/a;p,6'), fan(fcj) = xL f ™ ( {x r s, b - b',s,b), Xd 
Xd n{k) ( x P s > \b~ b'\), XdP = Xdp{x v s, \b- b'\). 



fan(fe) 



?! V 



*A A 




■in 



Figure 10. Total hadron-hadron eikonal, including enhanced diagram contributions. 



The first term on the r.h.s. in Fig. 10 is a simple Pomeron exchange 
(Xad)') tli e second graph contains a multi-Pomeron vertex which couples 
together any number (> 1) of "generalized fans" connected to the projectile 
hadron and any number (> 1) of those connected to the target; the third 
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graph subtracts Pomeron self-coupling contribution; the other terms on the 
r.h.s. in Fig. 10 correct for double counts of the same diagrams in the second 
graph. The latter can be verified using Figs. 7-9. 

Thus, hadron-hadron elastic scattering amplitude f a d(s, b) can be cal- 
culated using Eq. (3), with the simple Pomeron cikonal ^ ad {s,b) being 
replaced by X a °d^' k \ s ' b), corresponding to the full set of graphs of Fig. 10. 
To describe particle production one has to consider unitarity cuts of the 
amplitude f a d{s,b). Applying AGK cutting rules, 17 one can obtain a set 
of all cuts of diagrams of Fig. 10, with the eikonal x a °d cut ^' k ^ (s, b) ~ 
X a °d^' k \ s 'b), c in a similar form, i.e. as a contribution of one cut Pomeron 
exchange, the latter containing all relevant screening corrections (any num- 
ber of multi-Pomeron vertexes along the same cut Pomeron line, each ver- 
tex having any number but at least one uncut "fan" connected to the 
projectile or to the target); plus a graph with a multi-Pomeron vertex, 
which couples together any number (> 1) of cut "generalized fans" (includ- 
ing diffractive cuts) connected to the projectile hadron and any number 
(> 1) of those connected to the target, with all cut Pomeron lines being 
"dressed" with screening corrections in a similar way; minus double count- 
ing contributions. 10 This allows to obtain positively defined probabilities 
for various configurations of the interaction and to generate the latter via 
a MC method: starting from sampling (at a given impact parameter b 
and for given diffractive eigenstates j, k) a number of "elementary" parti- 
cle production processes - according to the eikonal X a °d" CUt ^^( s ^)' ano - 
choosing a particular sub-configuration for the latter - just one cut Pomeron 
exchange or a number of cut "generalized fans" emerging from the same 
multi-Pomeron vertex; in the latter case for each cut "fan" one generates 
further multi-Pomeron vertexes in an iterative fashion, according to the 
corresponding probabilities, until the process is terminated. Finally, like 
in the original linear scheme, 5 ' 6 ' 15 one performs energy-momentum sharing 
between constituent partons of the projectile and of the target, to which 
cut Pomerons are connected, and finishes with the treatment of particle 
production from all cut Pomerons. The described scheme easily generalizes 
to hadron- nucleus and nucleus- nucleus collisions, without introducing any 
additional parameters. There, screening effects arc enhanced by A-factors, 
as different Pomerons belonging to the same enhanced graph may couple 
to different nuclcons of the projectile (target). 



c Here the equality between the uncut and cut cikonals is approximate due to a somewhat 
different selection of neglected "chess-board" graphs in the two cases. 10 
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4. Cross Sections and Structure Functions 

We can also calculate screening corrections to sea quark and gluon PDFs 
at the initial virtuality scale Qq, which come only from "fan" -type dia- 
grams (neglecting Pomeron "loops"). Thus, parton momentum distribu- 
tions, being defined in the linear scheme by Eq. (8), are now described by 
the diagrams of Fig. 7, with the corresponding contributions (Eq. (13)) be- 
ing summed over hadron a diffractive eigenstates j (with weights C"A"), 
integrated over b', and with the down-most vertex being replaced by the 
Pomeron-parton coupling: 



Similarly one can express diffractive PDFs ff/^(x, xp, Qq) for a large 
rapidity gap process in deep inelastic scattering (with y gap = — Inxp) via 
the contribution of diffractive cuts of "fan" diagrams Xa°^ (s, Xp, b'). The 
latter is obtained cutting the diagrams of Fig. 7 according to the AGK rules 
and keeping contributions of cuts with the given rapidity gap, which yields 
again a recursive equation - Fig. 11: 



Correspondingly, the diffractive PDFs ff^(x, xp, Qq) are described by the 
diagrams of Fig. 11, with the corresponding contributions (Eq. (17)) being 
summed over hadron a diffractive eigenstates, integrated over b', and with 
the down-most vertex being replaced by the Pomeron-parton / coupling, 
i.e. with the eikonal Gp P being replaced by Gpf ft (cf. Eqs. (13), (16)). 
The obtained parton momentum distributions can be used to calculate 





(17) 
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Figure 11. Recursive equation for the diffractive cut \ a (s, xp, b') of "fan" diagrams, 
s = so/x. Cut Pomerons are marked by crosses. 

total and diffractive structure functions (SFs) F 2 , F^ 3 ^ as f 

F 2 {x,Q 2 ) = ejxf I/p (x,Q 2 ) (18) 
i=q,q 

F^ 3 \x r ,P,Q 2 ) = e]xfY»(x,x v ,Q% (19) 
i=q,q 

where (3 — x/xp and fi/ p (x,Q 2 ), ff^(x,xp,Q 2 ) are obtained evolving the 
input distributions fi/ p (x,Qo), fj'E(x,Xr,Q 2 ,) from Q 2 , to Q 2 (for valence 
quark PDFs q v (x, Q 2 ) a parameterized input (GRV94) 20 has been used). 

For the cutoff value Qq = 2 GeV 2 the Pomeron parameters have been 
finally fixed on the basis of experimental data on total hadron-hadron cross 
sections, elastic scattering slopes, and proton SFs F 2 , F^ 3 \ in particular, 
we had a P (0) = 0.205, a' P (0) = 0.09 GeV~ 2 , r 3P = 0.026 GeV" 1 , 7P = 0.5 
GeV" 1 . The obtained behavior of <t*£(s), F 2 (x 1 Q 2 ), and F 2 D(3) (x P , (3, Q 2 ) 
is shown in Figs. 12-14. In Figs. 12-13 we plot also a%£(s) and F 2 (x, Q 2 ) 
as calculated without enhanced diagram contributions, i.e. setting t"3p = 0. 

It is noteworthy that in the described scheme the contribution of semi- 
hard processes to the interaction eikonal can no longer be expressed in the 
usual factorized form, Eqs. (7), (9). Significant non-factorizable corrections 
come from graphs where at least one Pomeron is exchanged in parallel to the 
parton hard process, with the simplest example given by the 1st diagram 
on the r.h.s. in Fig. 3. In fact, such contributions are of importance to 
get a consistent description of hadron-hadron cross sections and hadron 
structure functions. At the same moment, due to the AGK cancellations 17 



'Strictly speaking, Eq. (19) is only appropriate in the small (3 limit; at finite the 
contribution of the so-called qq diffraction component becomes important. 19 
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3 4 5 

10 10 10 
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Figure 12. Total pp cross section calculated with and without enhanced diagram con- 
tributions - full and dashed curves correspondingly. The compilation of data is from 
Ref. 21. 




10 10 10 " 10 



Figure 13. Proton SF Fj(^,Q 2 ) calculated with and without enhanced graph correc- 
tions — full and dashed curves correspondingly. The data are from Refs. 22-24. 



the above-mentioned non-factorizable graphs do not contribute to inclusive 
high-p t jet spectra. Single inclusive particle cross sections are defined by 
diagrams of Fig. 15. In particular, parton jet spectra are thus obtained in 
the usual factorized form, being defined by the integrand in curly brackets in 
Eq. (9), with the corresponding PDFs ffu(x, Q 2 ) containing "fan" diagram 
screening corrections. 
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Figure 14. Diffractive proton SF F 2 (xp, /3,Q 2 ) in low /3 bins compared to prelimi- 
nary data of ZEUS collaboration (Ref. 25). 




Figure 15. Diagrams contributing to single inclusive cross sections; cut Pomerons are 
marked by crosses. V c {p c ) is the particle c emission vertex from a cut Pomeron. 



5. Results for Extensive Air Showers 

In general, various features of hadronic interactions contribute to air shower 
development in a rather non-trivial way. Nevertheless, basic EAS observ- 
ables are grossly defined by a limited number of macroscopic characteristics 
of hadron-air collisions: inelastic cross sections <7/£fLj, , so-called inelastic- 
ities Kf^_ ait (relative energy differences between the lab. energies of the 
initial and the most energetic final particles), and charged particle multi- 
plicities Nji^gfc. Indeed, the most fundamental EAS parameter, the position 
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of the shower maximum X max (atmospheric depth X, where the number of 
charged particles _/V e ± (X) reaches its maximal value), is mainly determined 
by C/Jlf'air and K 1 ™ 1 ^. In turn, the number of charged particles N e ±, mea- 
sured at a given observation level, is strongly correlated with X max . On 
the other hand, measured muon number iV M has a much smaller correlation 
with the shower maximum but depends significantly on A^_ air . 

In Figs. 16-18 the predictions of the new model for a^I 1 ^, K™^ l aiT , and 




Figure 16. Inelastic hadron-air cross sections calculated using QGSJET-II and QGS JET 
models - full and dashed curves correspondingly. 




° 4 io 4 io 6 io 8 io 10 

lab. energy, GeV 

Figure 17. Inelasticities of hadron-air interactions for QGSJET-II and QGSJET models 
- full and dashed curves correspondingly. 

^ht a ir are plotted in comparison with the results of the original QGSJET 
model. The observed different behavior results from a competition of two 
effects. First of all, non-linear screening corrections reduce the interaction 
cikonal, slow down the energy increase of hadron-air cross section, com- 
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Figure 18. Multiplicities of charged particles in hadron-air interactions for QGSJET-II 
and QGSJET models - full and dashed curves. 



pared to the linear scheme with the same Pomeron parameters, and sup- 
press secondary particle production. On the other hand, larger Pomeron 
intercept and steeper PDFs in QGSJET-II S lead to a faster energy increase 
of the latter quantities. At not too high energies the first effect dominates, 
resulting in smaller inelasticities and multiplicities. Nevertheless, in very 
high energy limit the influence of parton distributions prevails and the new 
model predicts larger values of K)f° l air and Nf^_ air . 

Sizably smaller inelasticities of the new model lead to a somewhat deeper 
position of the shower maximum, compared to the original QGSJET - 
Fig. 19. As the relative strength of non-linear effects is larger for nucleus- 
nucleus collisions, the difference in X max is stronger for nucleus- induced 
air showers. Although the effect is rather moderate at highest energies 
it changes drastically the interpretation of EAS data: while the predic- 
tions of the original QGSJET, being compared to, e.g., HIRES data, 26 are 
marginally consistent with the assumption of ultra high energy cosmic rays 
being only protons, this is no longer the case with the new model. 

The relative difference between QGSJET-II and QGSJET models for 
predicted muon number (E^ > 1 GeV) at sea level for vertical proton- 
and iron-induced showers is shown in Fig. 20. In the new model is 
significantly smaller, by as much as 30% at highest energies, due to the 
substantial reduction of N^_^ ir over a wide energy range (see Fig. 18). 

In general, applying the new model to EAS data reconstruction will 
change present conclusions concerning cosmic ray composition towards 



sin the original QGSJET model gluon PDFs are rather flat; xG(x, Q§) ~ x 1 -^ ) , with 
op(0) = 1.07, Ql = 4 GeV 2 . 
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Figure 19. Average position of the shower maximum for vertical proton- and iron- 
induced EAS as calculated with QGSJET-II and QGSJET models - full and dashed 
curves correspondingly, compared to the data of HIRES collaboration (points) - Ref. 26. 




10 10 10 

lab. energy (GeV) 



Figure 20. Relative JV M -difference (E^ > 1 GeV) at sea level between QGSJET-II and 
QGSJET models for vertical proton- and iron-induced EAS - full and dashed curves. 



heavier primaries. At the energies in the region of the "knee" of the CR 
spectrum (E ~ 10 15 — 10 16 eV) the obtained changes seem to go in the 
right direction to improve the agreement with experimental data. 27,28 How- 
ever, at highest energies the existing conflict between the composition re- 
sults obtained with fluorescence light-based measurements and with ground 
arrays 29 (much heavier composition in the latter case) will aggravate using 
the new model. Indeed, the first method is mainly based on the measured 
shower maximum position, whereas in the second case the composition re- 
sults arc obtained from studies of lateral muon densities at ground level. 
As the predicted reduction of JV M is much stronger compared to the corre- 
sponding effect for A max , the mismatch between the corresponding results 
is expected to increase. 
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6. Discussion 

The approach presented here offers a possibility to treat non-linear interac- 
tion effects explicitely in individual high energy hadron-hadron (hadron- 
nucleus, nucleus-nucleus) interactions. The latter are described phe- 
nomcnologically by means of enhanced (Pomeron-Pomeron interaction) di- 
agrams. A re-summation procedure has been worked out which allowed to 
take into account all essential enhanced graphs, both cut and uncut ones, 
to obtain positively defined probabilities for various configurations of in- 
elastic interactions, generally, of complicated topology, and to sample such 
configurations explicitely via a MC method. As an important feature of the 
proposed scheme, the contribution of semi-hard processes to the interaction 
eikonal contains an essential non-factorizable part. On the other hand, by 
virtue of the AGK cancellations, 17 the corresponding diagrams do not con- 
tribute to inclusive parton jet spectra. The latter are defined by the usual 
QCD factorization ansatz, with screening corrections contained in parton 
momentum distributions. 

Compared to other approaches, 7,30 the described scheme does not re- 
quire high parton densities to be reached; essential non-linear corrections 
to the interaction dynamics are consistently taken into account, providing 
a smooth transition between "dilute" (small energies, large impact param- 
eters) and "dense" regimes and approaching the saturation limit for "soft" 
(I? 2 1 < Qo) particle production in the latter case. On the other hand, 
the basic assumption of the scheme - "soft" process dominance of multi- 
Pomeron vertexes, poses restrictions on its application in the region of high 
parton densities. Indeed, after reaching parton density saturation at the 
chosen cutoff scale Q 2 ,, one has to account for corresponding effects in the 
perturbative (\q 2 \ > Qq) parton dynamics, i.e. to take into account "hard" 
Pomeron-Pomeron coupling. Such a treatment is missing here. 

Still, parton saturation is not observed in proton structure functions - 
Fig. 13. However, in our scheme the continuing increase of PDFs in the 
low x limit is partly due to the non-zero Pomeron slope, which leads to 
increasing contribution of proton periphery at small x. To investigate the 
problem in more detail we plot in Fig. 21 gluon momentum distribution 
for a given impact parameter x G(x, 6', Qq), the latter being defined by the 
integrand in curly brackets in Eq. (16). As seen from the Figure, at small 
b' the gluon density is indeed close to the saturation at the scale Qo m the 
small x limit, signalizing the need for perturbative non- linear corrections. 

Applying the approach to air shower simulations we obtained a siz- 
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Figure 21. Gluon PDF x G(x, b, Q 2 ) for a given impact parameter b (from up to bottom 
b = 0, 0.5, 1, 1.5 Fm) calculated with and without enhanced graph corrections — full and 
dashed curves correspondingly. 



able effect for the calculated EAS characteristics, with the shower maxi- 
mum getting deeper and with the EAS muon content being significantly re- 
duced. Accounting for perturbative screening corrections ("hard" Pomeron- 
Pomeron coupling) may enhance corresponding effects. On the other hand, 
the latter are only significant in the region of high parton densities, i.e. in 
the "black" region of small impact parameters. Thus, the likely outcome 
is a further reduction of the interaction multiplicity, without a significant 
effect on the cross sections and inelasticities, correspondingly, on the pre- 
dicted shower maximum position. Moreover, a possible loss of coherence in 
the leading hadron system, 31 which will only affect particle production in 
the fragmentation region, may have an opposite effect: increasing the in- 
elasticity but having the multiplicity essentially unchanged. Thus, the two 
effects will lead to further reduction of EAS muon number while having 
A max essentially unchanged, or even shifted upwards. Unfortunately, this 
will worsen the existing contradiction between fluorescence light-based and 
ground-based results on the cosmic ray composition. 29 There is a hope to 
clarify the situation with the data of Pierre Auger collaboration, 32 where 
both experimental techniques are available. The latter can provide an in- 
direct model consistency check at the highest CR energies. 
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